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We study the transport of charge due to polarons in a model of DNA which takes 
in account its 3D structure and the coupling of the electron wave function with 
the H-bond distortions and the twist motions of the base pairs. Perturbations 
of the ground states lead to moving polarons which travel long distances. The 
influence of parametric and structural disorder, due to the impact of the ambient, 
is considered, showing that the moving polarons survive to a certain degree of 
disorder. Comparison of the linear and tail analysis and the numerical results 
makes possible to obtain further information on the moving polaron properties. 



1. Introduction 

Charge transport along DNA is a subject of particular interest for two 
main reasons: on the one hand, it plays a fundamental role in biological 
functions as repair and biosynthesis; on the other hand, because of possible 
applications in molecular electronics and as molecular wires . Results on 
experimental DNA conductivity are controversial. It has been reported 
that it is good conductor 2 , insulator 3 and semiconductor 4 . A possible 
explanation, among others, for these striking differences can be the different 
long-range correlations of the DNA sequences 5 . 

In this work we propose a nonlinear mechanism for charge transport 
along DNA in the framework of the base pair picture ' , taking into account 
its spatial structure and the coupling of the spatial and electron variables. 
It turns out to be an efficient mechanism for charge transport which survives 
to a certain degree of diagonal and structural disorder. 
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2. Model 

We consider a semi-classical three-dimensional, tight-binding model for the 
DNA molecule, the sketch of the system being shown in Fig. 1. The lattice 




Figure 1. Sketch of the helical structure of the DNA model, the bases being represented 
by bullets. Geometrical parameters Ro, 9o, lo and the radial and angular variables r„ 
and respectively are indicated. 



oscillators are treated classically while the charge is described by a quantum 
system. The justification is that the nucleotides are large molecules with 
mass numbers of about 300. This also implies that molecular motions are 
small and slow compared to the one of the charge particle. The Hamiltonian 
of the whole system is given by H = H e \ + ff ra d + ff twist , where -ff ra d and 
_ff twist are, de facto, classical Hamiltonians and we can omit the hat on them. 
The Hamiltonians corresponding to the distances between nucleotides in 
each base pair ff r , and the twist motion -fftwist are given by 

ii - d ~z^w + 2 ' Ht ^-2^ ^7 + 2 ( j 

n n 

where {r n } represent the stretchings from the equilibrium distance, M the 
reduced mass or each base pair, f2 r the linear radial frequency, i} 
are the angles between two consecutive base pairs with respect to their 
equilibrium value 8 = 36°, J is the inertia moment, and fig, the twist 
linear frequency. 
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The electronic part is given by a tight-binding system of the form 
flei = E nH( n \ - V n -i, n \n - l)(n\ - V n+1 . n \n + l)(n| . (2) 

n 

In the state \n) the charge carrier is localized at the n th base pair. The 
quantities V n -i, n represent the nearest-neighbour transfer integrals along 
the base pairs and E n are the on-site matrix-elements. We write the elec- 
tronic state as \^S>) = ^2 n c n (t)\n), where c n (t) is the probability amplitude 
of finding the electron at the state |n). The interaction between the elec- 
tronic variables and the structure variables r n and 6 n . n -i arises from the 
dependence of the electronic parameters E n and V n ,n-\ on the spatial co- 
ordinates, given by E n = E n + kr n and V n , n -i = Vfo (1 — ad n , n -i), where 
d n ^ n -i is the first order Taylor expansion of the distances between con- 
secutive nucleotides, with respect to their equilibrium value. Is is given 

by 

R 

d n ,n-i - i— [(1 - cos(^o)) (r n + rn-!) + sin(6» ) R n ,n-i] • (3) 
k) 

Realistic parameters for DNA molecules are given in Refs. 8 > 9 ' 10 . We scale 
the time according to t — > O r t, which allow us to write the Hamiltonian 
in terms of dymensionless quantities 10 . The classical Hamiltonian is de- 
fined as -ffciass = {^\H\^), the dynamical equations for r n and # nj „_i are 
p n = Mf n = -dH c \ ass /dr n and p^„_ 1 = J0„,„-i = -S-Hciass/^n.n-i, 
while the evolution equations for the electron variables c„ are obtained 
from the Schrodinger equation iH(d 1 i'/dt) = H e i\SSf), which is equivalent to 
c„ = —{i/h)(dH c \ ass /dc n ). In this way we can obtain the scaled dynamical 
equations, which are 

irc„ = (E n + kr n )c n 

- (1 - ad n+ i. n ) c-n - (1 - ad n , n -i) c «-i (4) 
R 

r n = ~r„ - k |c„| 2 - a — (1 - cos6> ) 

x { [4+i c ™ + c„+i<] + [<c n _i + c„c*_!] } (5) 

R 2 

O n ,n-i = -^ 2 #„.„-i - aV sin 6»o [c* c„_i + CnC*^], (6) 

<o 

where r = ?if2 r /Vb represents the time scale separation between the fast 
electron motion and the slow bond vibrations. 

The values of the scaled parameters are 10 r = 0.2589, Q 2 = [0.709 - 
1.417] x 10~ 2 , V = 0.0823, R = 34.862 and l = 24.590, the time unit being 
~ 1.6 ps. Note that the time scales for the different variables differ in an 
order of magnitude: the fastest variables are the {c n }, with a characteristic 
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frequency of order 1/r ~ 4, followed by {r n } with the unity, and {# n ,n-i} 
with Q ~ 0.08 

There are no reliable values for the electron-radial and electron-twist 
coupling parameters, k and a. The criterion we have taken for their values 
is the consistency of the numerical simulations with the hypothesis, i.e., 
the deformations of the helix are small, and the linear approximations of 
the distance <i ni „_i and the trigonometric functions remain valid. Typical 
values are k ~ 1 and a ~ 0.002. 



3. Stationary states 

We suppose initially that r„ and 6 n ,n-\ are constant, i.e., the Born- 
Oppcnhcimcr approximation. This allows the obtention of expressions for 
them from Eqs. (5)-(6) which are inserted in Eq. (4), leading to a nonlinear 
Schrodingcr equation for the electronic amplitudes: 



\TC n 



El -k 2 \c n \ 2 -ka^ (1- cos o ) 



X { [ c *n+l c n + Cn+lC*] + [<C„_1 + C n C* n _^\ } ] C„ 

- ( 1 - ad n+ln ) c n+ i - ( 1 - ad„,„_i ) c„_i , (7) 

whith {g?„ ; „_i} given by Eq. (3), depending algebraically on the {c n } 
through {r„} and {0„,„_i}. 

To obtain stationary localized solutions we use a numerical method n > 12 . 
We substitute c n = <!>„ exp(— iEt/r) in Eq. (7), with $„ constants, and 
obtain a nonlinear difference system _E$ = A<j>, with $ = (<f>i, . . . , $jv)- 
Its solutions are attractors of the map: 

$ -> = i$/||i$|| , (8) 

|| • || being the quadratic norm. We start with a completely localized state 
at a site no, i.e., $„ = S n , no and apply the map above until <f>' = $. In this 
way we obtain stationary localized solutions and their energy E. 

In the ordered case, with E® = E 0} Vn, E can be made zero with the 
ansatz c„ — ► c n cxp(— lE^t/r). We obtain symmetric polarons with width 
of about 20 sites, combined with a local compression of order ~ 0.15 and 
a local unwinding of ~ 1° and energies E < — 2 ~ — 0.2cV. Considering 
static diagonal disorder with random E® the localization is enhanced, due 
to Anderson localization, the polaron being asymmetric, with a specific 
shape that depends on each particular disorder implementation. Similar 
results are obtained with structural disorder, i.e., random equilibrium val- 
ues i?" and 0°. The Floquet analysis shows that these polarons are linearly 
stable 10 . 
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4. Charge transport 

To activate the polaron motion we need to perturb the zero velocities of 
the ground state with a localized, spatially antisymmetric mode. This is 
obtained either with the time consuming pinning mode method 13 , cither 
with the simpler discrete gradient method, i.e., perturbing the variables 
r n with velocities parallel to {r n +i — r n -i} 14 . The main difference is 
that the first guarantees the mobility, while the second does not, but in 
practice works well very often. Values of the modulus of the kick velocity 
A r = (E 7 ^ )) 172 ~ °- 02 ' equivalent to a kinetic energy of ~ 200 meV are 
appropriate to obtain good mobility with low radiation. 

We obtain propagation of the localized electron amplitudes with con- 
stant velocity. A complementary compression and unwinding of the helix 
travels with it, while a localized oscillation of the angular variables re- 
mains at the initial positions. If we increase the diagonal disorder with 
E n G [—AE,AE], the charge transport still takes place, until values of 
AE ~ 0.05. Thereafter mobility becomes impossible and the polaron is 
pinned to the lattice. With disorder, however, the velocity of the first mo- 
mentum of the electronic occupation amplitude n c (t) = E„ n \c n (t)\ 2 is not 
uniform, as shown in Fig. 2. 

We also consider structural disorder with random distributions of the 
base pair spacings and angles around their ordered case values, with differ- 
ent mean standard deviations A. The conclusion is that, up to values of 
A = 10%, it does not significantly affect polaron mobility, although there 
is a slight reduction of its velocity. 

It is interesting to remark that for the values of a coherent with our small 
amplitude hypotheses it is not possible to move the polaron by kicking the 
angular variables. These results are coherent with the fact that the Floquet 
analysis shows that the pinning mode appears only in the radial variables 
and its eigenvalue separates from the optimal value around 1 when the 
disorder is increased 10 . 



5. Linear and tail analysis 

Although our system is studied in the nonlinear regime, and the results 
obtained throughout this work are essentially nonlinear, there is a number 
of linear techniques that can be applied as a reference and as a tool for 
obtaining useful information. They are the study of the linear system itself 
and the tail analysis. 
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Figure 2. Diagonal disorder. The position of the center of the electron breather as a 
function of time and for different amounts of disorder. Full line (Ordered case), short 
dashed line (AE = 0.025), dashed-dotted line (AE = 0.050), dotted line (AE = 0.100) 
and long dashed line (AE = 0.500). 

5.1. The linear system 

If in Eqs. (4)-(6) we cancel out the nonlinear terms, we obtain: 

i T C n — C n -\-i C n — 1 , T n = r n , Qnn—l = ^ @nn—l • (9) 

In the first equation the terms E n , which are all equal to some value E 
in an homogeneous system, have been eliminated through the ansatz c„ — > 
c n exp(— iEo t/r). The system becomes decoupled and the variables r n 
and Q nn -\ correspond to independent linear oscillators with frequencies 
w r = 1 and £1 = 0.0842, respectively, in the scaled units. To obtain the 
ground state we substitute in Eq. (9) c n (t) = (j) n exp(— iEt/r), 4> n being 
time independent. We obtain the stationary discrete Schrodinger equation: 
E(j)n — — </>n+i — 4>n-i- Substitution of the linear modes <p n — exp(iqn) 
leads to: E = —2 cosq. Therefore, the linear energy spectrum runs from 
—2 to 2. The minimum energy, E = — 2, corresponds to the wave vector 
q = 0, i.e., all the oscillators vibrating in phase. The nonlinear ground 
states described above have E < —2 and derive from this mode. 
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5.2. Tail analysis 

For sites far enough from the polaron center, which can be only a few 
sites, we can still apply the linear analysis. Now, the substitution of the 
tail mode c n (t) = </) n exp(-( n — iqn), with n > and £ > 0, leads to: 
E = — 2 cos(g) cosh(£) and 2 sin(g) sinh(£) = 0. The second equation, 
implies that only two wave vectors are possible, q = 0, with negative energy, 
and q = it, with positive one. The energy of the first, which is the one 
that appears throughout this paper because it is movable, is given by: 
E = — 2cos(£) < —2. Therefore, the nonlinearity, produces localized states 
with larger values of the frequency w p — \E\/t outside the linear spectrum. 
The distance of the energy calculated in the full system from E q= o = — 2, 
is a measure of the degree of nonlinearity. It is interesting to relate the 
values of the energies with the polaron breadth. Here, we define it loosely 
as three times the number of sites necessary for |^ n | 2 to decrease to a 5%. 
The factor 3 is there to allow for both and for the nonlinear center of the 
polaron. This leads to An ~ 4.5/£. An approximate table is: 





2 


1 


0.5 


0.15 


0.075 


An 


2.25 


4.5 


9 


30 


60 


E 


-7 


-3 


-2.25 


-2.02 


-2.006 



where the dimensionless E units are equivalent to 0.1 eV. This makeshift 
method gives a good estimate of the breadth of the excitation and produces 
a good fitting a few units far from its center. The numerically calculated 
polarons found in this work correspond to An ~ 20. 

The tail analysis can also be applied to the moving polaron, far enough 
from its center. Let us propose a localized traveling wave of the form: 
c„ = exp{— £ (n + v t) — i (qn + w p t)}, which represents a traveling localized 
wave moving to the left behind the polaron, with positive £, q, and velocity 
v. Substitution in Eq. (9) leads to the following equations: E = w p t = 
—2 cos(g) cosh(£) and t£v = 2 sin(q) sinh(£). They are undetermined, as 
there are four unknowns: E, q, £ and v, and do not allow us to determine 
the velocity. This is little surprise, as the polaron velocity depends on the 
energy given with the perturbation to the static one, and this is done to 
the spatial variables. However, these equations can be used to obtain an 
estimate of the moving polaron energy and wave number, using the decay 
length £ corresponding to the static one, as there is no appreciable change 
of shape, and the velocity observed within the simulations. The conclusion 
is that there is only a negligible increase of the energy, i.e., the energy 
given by the kick is stored in the spatial variables. The fact that there is 
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a imaginary part of the energy —i^vr is, as it is known, consequence of 
the fact that we are dealing with only a part of the system 15 , in which the 
energy is decreasing as the polaron moves to the left. 

6. Conclusions 

We have considered a model for charge transport along DNA, taking into 
account the 3D structure. We obtain linearly stable, localized stationary 
states. We are able to move them using the pinning mode and the dis- 
crete gradient methods. For the ordered case, the translational velocity 
is constant, while for the disordered case the electron motion is irregular 
and if the disorder is high enough it is impossible. We have performed a 
linear and tail analysis of the system, which makes possible to obtain the 
phonon spectrum, measure the degree of nonlinearity, obtain estimates of 
the polaron breath and energy. It shows that the charge transports very 
little energy, being most of it associate with the spatial coordinates. The 
main conclusion of the whole work can be that our proposed mechanism is 
an efficient means for charge transport along ordered and disordered DNA. 
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